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Abstract The analysis of non-radiative sources of static or time-dependent gravita- 
tional fields in the Solar System is crucial to accurately estimate the free-fall orbits of 
the LISA space mission. In particular, we take into account the gravitational effects of 
Interplanetary Dust (ID) on the spacecraft trajectories. The perturbing gravitational 
field has been calculated for some ID density distributions that fit the observed zodia- 
cal light. Then we integrated the Gauss planetary equations to get the deviations from 
the LISA keplerian orbits around the Sun. This analysis can be eventually extended 



M. Cerdonio and F. Marzari 

Department of Physics, University of Padova and INFN Padova, via Marzolo 8, 1-35131 Padova, 
Italy 

F. De Marchi 

Department of Physics, University of Trento and INFN Trento, via Sommarive 14, 1-38123 
Povo (Trento), Italy 

R. De Pietri 

Department of Physics, University of Parma and INFN Parma, 1-43100 Parma, Italy 
Ph. Jetzer 

Institute of Theoretical Physics, University of Zurich, Winterthurerstrasse 190, CH-8057 
Zurich, Switzerland 

G. Mazzolo* 

Max Planck Institut fur Gravitationsphysik, Callinstrasse 38, 30167 Hannover, Germany 
*To whom address correspondence: giulio.mazzolo@aei.mpg.de 

A. Ortolan 

INFN Laboratori Nazionali di Lcgnaro, Vialc dell'Universita 2, 1-35020 Legnaro (Padova), 
Italy 

M. Sereno 

Institute of Theoretical Physics, University of Zurich, Winterthurerstrasse 190, CH-8057 
Zurich, Switzerland and Department of Physics, Politccnico di Torino, Corso Duca dcgli 
Abruzzi 24, 10129 Torino, Italy 



2 



to Local Dark Matter (LDM), as gravitational fields are expected to be similar for ID 
and LDM distributions. Under some strong assumptions on the displacement noise at 
very low frequency, the Doppler data collected during the whole LISA mission could 
provide upper limits on ID and LDM densities. 

Keywords LISA ■ interplanetary dust • dark matter • gravitational waves 



1 Introduction 

LISA (Laser Interferometer Space Antenna) is a joint space mission by ESA and NASA 
which is planned to be launched at the end of the next decade. It consists of three 
identical free-falling satellites, orbiting around the Sun and marking the vertices of a 
nearly equilateral triangle with ~ 5 • 10 6 km (1/30 AU) sides, located 20° behind the 
Earth and lying in a plane that makes an angle of 60° with the ecliptic |2j. LISA target 
is the detection of gravitational waves (GWs) through the measure of the relative and 
differential motions between the spacecrafts. Among the astrophysical goals of LISA is 
to detect GWs originated by events like black holes coalescence or capture of compact 
objects by black holes [3]. This requires a strain sensitivity of 10~ 20 Hz -1 / 2 < S^ 2 < 
1CT 17 Hz~ 1/2 in the 1CT 4 to 1CT 1 Hz frequency band range. 

However, LISA reference masses interact also with time-dependent (e.g. planets 
and their satellites, quadrupolar pulsation of the Sun, etc.) and static (e.g. ID and 
LDM) local gravitational fields, and therefore they depart from ideal unperturbed 
orbits around the Sun. In this paper, we focused on the perturbing effects caused by 
the static components ID and LDM. As the discriminating feature between ID and 
LDM concerns only their coupling with the electromagnetic field, we expect that ID 
and LDM will produce similar gravitational forces on reference masses^] In fact, only 
ID reflects the solar light and can be studied through observations of zodiacal light [5] , 
while LDM is "dark" and its presence can be investigated by means of gravitational 
perturbations induced on orbiting bodies. It is worth mentioning that other possibilities 
for detection of DM are the searches for particular electro-weak processes beyond the 
Standard Model at particle colliders. Thus we assume that ID and LDM induce identical 
gravitational effects on LISA orbits, with the irrelevant difference that the first is 
luminous and the second is dark. 

As we will show, the linear perturbation theory holds, being other gravitational 
forces acting on LISA reference masses much smaller than the Sun pull; under this 
approximation, small perturbations to ideal keplerian motions induced by different 
sources in the Solar System can be studied separately. 



2 Interplanetary dust 

The ID cloud is a dust composed by grains with typical sizes of 10 — 100 /im pervading 
the Solar System. The distribution of ID particles is studied by the observations of 
solar light scattering (i.e. the zodiacal light confined to the ecliptic plane) and by their 

1 This is correct only in Newtonian physics: in fact, DM and ID are expected to be gravi- 
tationally bound to the Galaxy and to the Solar System respectively; therefore, while gravi- 
toclcctric effects arc identical for the two, that would not be the case for the gravitomagnctic 
ones, i.e., effects depending on the speed of the considered objects. 
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thermal emission, which is the dominant component of night-sky light in the 5 — 50 /im 
wavelength Q3]. Physical models of ID density distribution p should account for its 
symmetries, that can be easily described in the Solar System Baricentric (SSB) refer- 
ence frame (x,y, z): i) invariance under rotation about the z axis, and ii) invariance 
under reflections in the (x,y) ecliptic plane [^] Usually one also assumes a static dis- 
tribution and a power-law radial profile. As a consequence, all the proposed models 
factorize p into two functions [5] 

p(r,/3 ) = p o (^) Q /(/?0), (1) 

where r is the distance from the Sun, (3q is the helioecliptic latitude, pQ ~ 9.6 x 1CP 20 
kg/m 3 is the density value at tq = 1 AU (Earth's orbit), as estimated by the integration 
of the meteoroid mass distribution far away from the Earth, in order to avoid its 
gravitational attraction on ID particles, and / is a given function. The typical value of 
the parameter a, determined from zodiacal light photometer measurements on Helios 
1 and 2 [10], is 1.3. 

The expression of / (/?©) is still uncertain; however, the analytical formula which 
best reproduces the observations of zodiacal light is the so-called ellipsoid model 

p(r,fe) = 1 (2) 

[1 + (7,gsm/?0) j 

where je = V a? — 6 2 /b, a and b are the semi-major and semi-minor axes of an oblate 
ellipsoid respectively. Beyond ~ 3 AU in the ecliptic plane and ~ 1.5 AU off the ecliptic 
plane, no reliable density values can be obtained from the zodiacal light. Therefore we 
assume a = 3 AU and b = 1.5 AU for the semi-axes and 7_e = y/3 [5]. It is worth 
noticing that the total ID mass amount inside the considered region is of the order of 
10 17 kg, i.e., « 10~ 8 Af© . 

The ellipsoid model of ID density depends on two parameters, a and 7g, and, in 
order to numerically study gravitational effects of ID on LISA orbits, we will consider 
four cases: i) spherical homogeneous distribution (a = 0,7s = 0), ii) spherical distri- 
bution with power-law density profile (a = 1.3, "/e — 0), hi) ellipsoidal homogeneous 
distribution (a = 0,7^; = v3)j an d iv) ellipsoidal distribution with power-law density 
profile (q = 1.3,7b = V3). 



3 Method 

As the first step we define the unperturbed LISA orbit^] The cartesian coordinates 
of each LISA reference mass (k = 1, 2, 3) are related to the keplerian orbital elements 
defined in the SSB through the following equations: 



2 The distribution is not symmetric with respect to the ecliptic plane, but to one defined by 
the total angular momentum of the Solar System. The differences between the two planes are 
negligible for our scopes. 

3 In this paper we neglect post-Newtonian and relativistic effects. 
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x k{~t) — a k [(cos J7fc cosa;^ — sin sinajj, cosi^,) (cos V'fc(i) — e fc)] 
+a& (— cos J7fc sintjj, — sin coscj/. cosi^) *Jl — e? sin ■j/'fe (t) 

/ J/fe(£) — a fc [(cos cosajj; + cos sina;^ cosife) (cosipk(t) — e^)] 
+a& (— sin svau^ + cos cosloj- cosi/.) -\/l — e 2 , sin ip^ (t) 



(3) 



sin Wfc sin ij. (cos tp^ (t) — e&) + cos ojj, sin ij. v/l"— e? sin ^ (i) 



where aj. is the semi-major axis, e/. the eccentricity, ij, the inclination of the orbit 
with respect to the ecliptic plane, to^ the argument of periapsis, i?^ the longitude of 
ascending node and ipkW) the eccentric anomaly [BJ- We chose the initial conditions 
providing rigid triangular configuration of side I ~ 5 x 10 6 km ~ 1/30 AU 



' a k = 1 AU 

e fe=( 1 + ^+l^) 1/2 - 1 

tanii, = ri r- 

, k l+?/V3 ( 4 ) 

. , TT 

^k — 2 

% = §7r(fe-l)-f 

t M k (t) = 2-nt - §7r(fc - 1) - 7T, 

where £ = i/(2a/») ~ 1.67 x 10~ 2 and M&(t) is the mean anomaly. Such orbits minimize 
the variation of the inter-spacecraft distance between the i-th and j'-th satellite AL^ j 
= \vi(t) — r j- (i ) | |12| . which turns out to be independent of time at the first order in £ 
[8]. This requirement must be fulfilled in order to operate LISA successfully. In fact, the 
measured quantities for the GW searches are the differential relative motion between 
the couples and of the three reference masses, &L^ j i = AL^ j — ALj i 

(see Fig. {l}). Each couple can be regarded as an unequal arm interferometer and the 
doppler shift induced by the relative differential motions can be measured by a suitable 
Time Delay Interferometry (TDI) [T] with a strain sensitivity in Fig. |2j) [?]. 

The keplerian motions of the LISA reference masses are periodic. We can distinguish 
their harmonics by looking at the modulus of the Fourier transform of 5L, j,;(t) I as 
in Fig. (j3|, where we have considered a finite observation time T and, to overcome the 
spectral leakage, we have applied a Blackman tapering function [15] . 

It is worth noticing the absence of the harmonics at 3, 6, 9 . . . y _1 due to the 
presence of a discrete symmetry in the equations of motion of the constellation. In 
fact, as a consequence of the initial conditions in Eq. Q, the LISA triangle rotates 
as a "quasi-rigid body" with period of one year and, after integer multiples of one 
third of year, the dynamical configuration is identical to the initial one under any 

4 The Fourier transform of F(t) is as usual F(u) = J + °° F(t) exp(—iuit)dt. 
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Fig. 1 Left, relative motion of 1 and 2 reference masses AL\^ as a function of time. Right, 
differential relative motion <5-Li,2,3 between couples 1-2 and 2-3 of reference masses as a function 
of time. For different couples, plots are identical but shifted in time. 




cyclic permutation of the reference mass indices. However, such a feature is only of 
theoretical interest, because slightly different initial conditions, due for instance to the 
unavoidable injection errors of LISA spacecrafts in their orbits (position ~ 2 Km and 
velocity ~ 2 mm/s, [17]), will produce harmonics of the 3 y _1 frequency. 

To study the effects of ID on LISA orbits, we apply the perturbation theory in 
the six-dimensional space of parameters and use the Gauss planetary equations. Such 
equations provide time evolution of the orbital parameters under a generic perturbing 
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Fig. 3 |<5Li, 2 ,3(o;)| windowed by a Blackman function with T = 30 y (1 y~ x « 3.1709792 x 
10" 8 Hz). 



acceleration field 7(r) and read 



( da k 

dt 



dt.-,. 



dt n 



n k v/l— e 



[ efe A r sin0 fc + A t (&)] 



^ = ^== A n (^)cos( Wfc + tf fc ) 

n k a k ^l-e 2 k 



d.Q k 
dt 



du) k 
dt 



dM k 
dt 



- sin l A M£)sinK + fl fe ) 
a fc sin i fc ^/ 1 — e fc 



a-e? 

- cos it + - 31 5 

K dt n k a k e k 



-A r COS 0fc + A, 



sm 



71 fe (-^) - J\ 

K ri k a k \a k ) V 



2 dtjfc 



(5) 



where = a^, (1 — e|) is the semi-latus rectum, 0j. is the true anomaly, nj. = */ 47r 2 /a^ 
is the keplerian mean motion and A r ,At and A n are the components of 7 along the 
versors f — r/ ||r|| in the radial direction, t, orthogonal to f in the osculating plane and 
in the direction of r, and h = f X t [6]. 

For the generic point r = (x, y, z) close to the ecliptic plane (ij. ~ 10~ 2 rad Vfc), 
the accelerating field -y(r) has been calculated from the internal gravitational potential 
of the distribution in Eq. Q, 
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OO 



R z 



r + a 



l 



R 2 T + b 2 



(t + a 2 ) Vt + b 2 



-dr. 



(6) 



where R = ^ x 2 + y 2 , ip = 2np Gr$a 3 ~ a yjl - y|/ (l + J%) / (2 - a) and G is the 
Newton gravitational constant. The gravitational potentials of the considered distribu- 
tions can be easily obtained for the appropriate values of a and je- 

We then numerically integrated the Gauss equations for the three LISA satellites 
adopting different ID distribution;^] In order to check the effectiveness of the linear 
perturbation theory of the studied perturbations, we have compared our results with 
the numerical solutions obtained by keeping the orbital elements on the RHS unper- 
turbed, and they matched exactly within the integration time (30 years). In addition, 
the accuracy of the numerical integrations (~ 10~ 30 ) have been verified by comparison 
with analytical solutions of Gauss equations (approximated to the fourth order in e 
and with unperturbed orbital parameters on the RHS) for a spherical homogeneous 
distribution of matter [7] and the agreement has been satisfactory. To get the per- 
turbed orbits, and, therefore, the perturbative effects of ID on LISA constellation, we 
substituted the unperturbed orbital parameters appearing in Eq. (|3| with the solutions 
of Eq. {5}. 



4 Results 

In Fig. Q, we report the plot representing 

«W=«^s-«X ( 7 ) 

i.e., the difference between perturbed and unperturbed differential motions of the con- 
stellation for the four considered distributions. 

We note that the perturbative effects of ID are very small, in fact the LISA constel- 
lation opens with an amplitude of the order of 10 2 fxm after 5 complete orbits. Fig. |5| 
shows the modulus of the Fourier transform of the difference between perturbed and 
unperturbed differential motions 5/1,2,3, f° r the ellipsoidal ID distribution with radial 
profile, and we see that ID enhances resonance peaks, in particular the 2 y _1 harmonic. 
In addition, new harmonics corresponding to integer multiples of 3 y _1 appear, as a 
consequence the ID gravitational field that breaks the permutation symmetry of the 
equations of motion. 

Therefore, our analysis for the different ID distributions listed at the end of Sec. 
Q shows that curves in Fig. Q and Fig. Q change their amplitudes of a factor w 2.5 
by varying 7^; from \/3 to 0, while they do not depend significantly on a. Additionally, 
it is easy to show that LISA sensitivity curve is not significantly affected by the ID 
perturbing effects in the I0~ 4 — I0 _1 Hz frequency band. 

5 We made use of the Mathematica 7 numerical integrator NDSolve; see the documentation 
webpage |http://reference.wolfram.com/mathematica/r ef/NDSol ve.html| 
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Fig. 4 (color on line) Time evolution of (5/1,2,3, for the spherical (je = 0) distribution (thin 
dashed line, a = 0, thin solid line, a = 1.3) and ellipsoidal (7^; = V3) (thick dashed line, 
a = 0, thin solid line, a = 1.3) distribution. It is worth noticing that ID induces effects of the 
order of 10 2 fim, 10 12 times smaller than those due to the keplerian differential motions. 
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Fig. 5 Difference between <5£f e ^"*3 and <5L" n*^*" spectra of frequencies after 30 orbits, win- 
dowed with Blackman function, obtained with a = 1.3 and 7^; = v3; we note the presence of 
the harmonics corresponding to the integer multiples of 3 y~ 1 , contrary to Fig |3|. 
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4.1 Optimal ID signal resolution 

It is of some interest to investigate the problem of resolving the contributions to the 
differential relative motions due to Sun and ID by means of optimal filtering. From 
the theory of signals resolution we known that if a noisy signal v(t) can be a linear 
combination of two given signals f(t) and g(t) of unknown amplitudes A and B 

v(t) = Af(t) + Bg(t) + n(t), (8) 

where n(t) is a zero mean gaussian stochastic process with correlation (n(t)n(i')} — 
SoS(t — t'), the four following hypothesis are possible: 

1. Hq: neither signal A nor signal B is present; 

2. H\: signal A alone is present; 

3. H?: signal B alone is present; 

4. H%: both signal A and B are present, 

which are to be verified by comparing the expectation values of A and B, with zero 
(test of hypothesis). The best estimate of A and B turns out to be |11| 

B = ^h? I W)-Xf{t))v{t)dt, (9) 



with standard deviations 



S, 



1/2 







-y/2 (1 — A 2 ) 



(10) 



where A = L f(t)g(t)dt, Sq is the noise power density and f(t) and g(t) are normalized 

rT 9 rT 9 

to unit energy, J Q / (t)dt = J Q g (t)dt = 1. 

In our case, v(t) represents the perturbed differential motions, Af(t) and Bg(t) are 
the contributes due to Sun and ID respectively. We stress that f(t), g(t), A and B 
are almost independent of the initial conditions of LISA reference masses. It is worth 
noticing that the two contributions to the perturbed motion are almost orthogonal, 
(A ~ 10~ 2 ) as can be seen in Fig. (JsJ) , where the normalized functions f(t) and g(t) 
are plotted. 

As a crude estimation, a strain sensitivity ~ 10~ n 1/VSz at the 2 y _1 (~ 
6.4 • 10 -8 Hz) frequency (the most powerful ID harmonic), is required to achieve a 
unitary signal-to-noise ratio (SNR = A/a^ = B/a^.) in a 5 years observation time. 
However, in the literature there exist extrapolations of the strain sensitivity Sh only 
up to 10 -6 Hz [3] and the measure of such a small signal amplitude depends on the 
noise level of LISA at very low frequencies (10~ 7 Hz), which is still largely unknown. 



5 Discussion 



We have calculated the ID perturbing effects on the LISA differential motions. Using 
the estimated ID density, we found a continuous opening of the LISA orbits of the 
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Fig. 6 (color on line) Time evolution of the two normalized functions f(t) (blue solid line) 
and g(t) (red dashed line) for the case of an ellipsoidal distribution with radial profile. 



order of 10 2 /im after 5 years at 2 y _1 frequency. Moreover, due to the particular LISA 
orbits, which are nearly circular and very close to the ecliptic plane, we found similar 
effects for the distributions we have taken into account. In fact, the ID tidal forces vary 
significantly over a scale much larger than the LISA triangle sides. 

All the results we presented for ID also hold in the presence of LDM. In fact, we 
have shown that the numerical solutions of Eq. (J5|, for the LISA orbits, are almost 
independent on a and je parameters within the investigated range and of the injection 
errors. In addition, the perturbations on LISA relative differential motions SLijk scale 
linearly with po- Thus we can account for LDM just by a rescale of po: for instance, 
if we assume that the LDM density value is close to the average galactic dark matter 
(GDM) density value, po,GDA/ — 5- 1CP 22 kg/m 3 as obtained from the galaxy rotation 
curve, the effect of DM is expected to be 0.5% that of ID. 

As a consequence, LISA could provide interesting upper limits on both Pojd and 
Po,LDM, depending on low frequencies noise and po value, by means of direct gravita- 
tional field measurements instead of observations of zodiacal light in the case of ID. At 
present time, the best upper limits po,LDM < 3 ■ 1CP 16 kg/m 3 are based on the study 
of the precession of the perihelions of Mercury, Earth and Mars |13| . 

However, we stress that a thorough study of LISA displacement noise below 1 pHz, 
including local gravitational field fluctuations, thruster noise, orbit determination and 
injection errors [5] etc., is required to establish the relevance of the ID effects in LISA 
physics. 



6 Conclusions 

The study of the Solar System gravitational field acting on LISA reference masses is of 
some relevance also for the detection of GWs. In fact, the non-radiative gravitational 
perturbations on LISA keplerian motions must be subtracted at the highest accuracy 
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from the relative differential motion in order to measure the GW contributions. In this 
paper we have calculated the time evolutions of the orbital parameters due to some ID 
distributions and estimated the opening induced on the LISA constellation, « I0 2 /jm 
after 5 years. The Fourier transform of the perturbation of the differential relative 
motion has shown an enhancement of the resonances characterizing the unperturbed 
spectra, in particular the peak at 2 y -1 frequency. As a consequence of the LISA orbits, 
such effects are similar for the studied distributions of ID and do not affect the LISA 
sensitivity band for GW detection. 

On the other hand, the discrimination of very small perturbations from the keple- 
rian differential relative motion depends crucially on the low frequency displacement 
noise, that in the ideal case should be ~ 10~ n l/\/Hz at 2 y _1 frequency to detect 
ID at the density ~ 10~ 20 kg/m measured by the zodiacal light. Unfortunately, a 
reliable estimate of this noise level at very low frequency is still lacking and requires 
futher investigations. 

Finally, we investigate the possibility of constraining the LDM density with LISA. 
Even though gravitational perturbations due to ID and LDM are undistinguishable, 
the study of the deviations from the keplerian orbits of LISA reference masses could 
provide interesting upper limits on LDM density. 
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